Bounded Gain of Energy on the Breathing Circle Billiard 
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Abstract 

> ' 

■ The Breathing Circle is a 2-dimensional generalization of the Fermi Accelerator. It is shown that the 

| billiard map associated to this model has invariant curves in phase space, implying that any particle 

. will have bounded gain of energy. 

H 

00 . 

Os ■ 1 Introduction 

la: 

The breathing circle is the time-dependent plane region bounded, instantaneously, by a circle whose 
q ■ radius is a periodic function of time. The billiard problem on the breathing circle consists on the free 

motion of a point particle inside this region, colliding elastically with the moving boundary. 

As J. Koiller, R. Markarian and ourselves have showed in ||, time varying billiard computations can 
be performed in the same way as done on rigid billiards. But, beside the usual coordinates describing 
static billiards (for the impact point on the boundary and the direction of the movement), one must also 
introduce the energy and time. If, for each fixed time, the frozen boundary is a strictly convex plane 
curve, the associated time-dependent billiards can be modeled by a 4-dimensional, volume preserving 
diffeomorphism which maps successive impacts with the moving boundary. 

In the special case of the breathing circle, the angular momentum is conserved, implying a reduction 
of the model to a 2-dimensional diffeomorphism. Then, using a corollary of Herman's Theoreme des 
Courbes Translatees Q we prove that if the motion of the moving boundary is sufficiently smooth then, 
for any admissible initial condition, the particle will move with bounded velocity and therefore will have 
a bounded gain of energy. 

A similar case was studied by Levi in §1), where he shows examples of pulsating soft circular billiards 
whose energies stay bounded for all time. 

The study of time-dependent billiards is motivated by two main questions. From one side, it is a natural 
extension to higher dimensions of the Fermi Accelerator, studied for instance in 0, M and ||. From 
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another side, among many branches of physics, it is of special interest on the study of the motion of 
particles inside collectively excited nuclei. On the attempts to understand the origin of dissipation of 
collective motion in finite Fermi systems (see, for instance, and models were constructed where 
the classical part corresponds to time-dependent billiards, especially with multipole resonances (with the 
breathing circle corresponding to the monopole mode). 



2 The breathing circle billiard map 

Let R(t) be a strictly positive T-periodic C k function, k > 2 and T(t) = {x 2 + y 2 = R 2 {t)} be the 
(moving) boundary of the breathing circle. Suppose that at an instant t n a particle is on T(t n ) at a point 
(x n ,y n ). Let the unitary tangent vector at that point be f„ and the unitary outward normal vector be 
fj n . The point on T{t n ) is then given by (x n ,y n ) = R n = R n f] n = R(t n )fj n . Suppose that it leaves the 
circle from (x n ,y n ) with velocity v n . 

It travels on a straight line and hits the boundary again at an instant t n+ i at {x n+ i,y n+ {) = R n +i with 

Rn+l — Rn + (t n +l — tn)V n (1) 

and will then rebound with velocity 

dR(t) 



V n +1 — V n - 2 < V n ,Vn+l > f] n+1 + 2 



(It 



Vn+1 ■ (2) 



t=t, 



+i 



Let L — v x R denote the angular momentum with respect to the center of the breathing circle, then 

L n +1 = Vn+l X Rn+l = V n X R n +1 = V n X R n = L n . 

This shows that the angular momentum is conserved at impacts with the moving boundary. As it is also 
obviously conserved between impacts (since the particles moves on a straight line with constant velocity), 
angular momentum is conserved along all the movement. A first easy consequence of this conservation 
is that a particle never changes orientation with respect to the breathing circle, i.e., it always rotates in 
the same direction. In other words, if we introduce the notation L n = R n v n ■ f n , is is easy to verify that 
L n is conserved. 

Another consequence of the conservation of angular momentum is that a particle never stops if its angular 
momentum is different from zero. This is because if v n = \v n \, as L n = Lq then R n v n > \R n v n • f„| 
/.:: and v n > 7^ > where R= max {R(t),t G [0,T)} . 

It also follows that the traveling time of a particle (between two impacts with the moving boundary) is 
bounded, because 

_ distance 2R 2R 2 

tn+l — t n — — < < — — (3) 

velocity v n Lq 

for every n. 

We introduce the variable I n = —R n v n ■ i) n . It is, in a certain sense, a natural variable, because it is the 
normal counterpart of the conserved quantity L n . 
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For given Lq ^ and Iq, to initial impact variables, let I\ and t\ be the values of / and t at the next 

^ 2 (t)| t=t „ 



impact with the breathing circle. Let w n = 4?R 2 (t)\ t 



Lemma 1 For any fixed Lq ^ and for every (Jo, to) with Iq + ^ > 0, the map M : (Io,t Q ) i— > (Ix,tx) 
is defined by the formulae: 

I 1 =I Q - Wl+ (^i^l(t 1 -t )-2I Q ^ (4) 
(*i ~ *o) ( X °^2 /0 (*i ~ *o) ~ 2J o) = iff ~ #0 (5) 



Proof: Equation (||) above is easily obtained after computing |i?i| 2 from (Q). Now, if the particle moves 
inside the breathing circle region, it's normal velocity at the impacts must be bigger than the normal 
velocity of the boundary, i.e., —Vq ■ fjo + R'(to) > or, using the variables introduced above, Iq + ^ > 0. 
This condition implies that (^|) has a solution t\ > to and thus the map is well defined under this 
assumption. (In || we have in fact proved a little more: Iq + ^ > implies ii + ^ > 0.) 

On the other hand, from (^) we have that (t\ — to)vQ 1 = {R\ — Ro) ' ^o- This relation together with 
the rebound formula (||) yields ([|). ■ 

Formulae (^) and (||) above are invariant under the translations t + nT, where T is the period of R(t) 
and n is any integer. So, we can take t (mod T) £ S 1 and then the domain of M. is contained in the 
cylinder JJxS 1 . It is worthwhile to point out that, although time-dependent billiards are described by 
four variables, the breathing circle only takes two, namely / and t (mod T), the other two independent 
variables being the (conserved) angular momentum and its conjugate variable. Moreover, as Lq appears 
only squared in the calculations, we can suppose Lq > 0. The case Lq = corresponds to the movement 
on the diameter of the breathing circle and is equivalent to the 1-dimensional Fermi accelerator, studied, 
for instance, in Q, and |j. 

Now let (It, fa) — M(Io,t ) and define (7_i,t_i) by M(I-t,t-i) = (Io,t ). Let 

A x = {{IqM I *o + y > and/i + y =0} 
A-x = {(J ,*o) \Io + Y >0 and/_ 1 + ^i=0}. 
(With our assumptions on R(t), the Lebesgue measure on M 2 of these sets is zero.) 
Let D = R x S 1 \ (AxUA-x)- 

Lemma 2 If R(t) is a strictly positive T -periodic C k function, k > 2 then for each fixed Lq ^ the map 
M. : D — > D, given by and (j^j, is a C k ~ 2 diffeomorphism, preserving the area d/i = ^ 2 dldt. 



Proof: the proof follows from Theorem 1 in 
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3 Changing coordinates: Some properties of the billiard map 



From the section above, it seems natural to introduce the variable J = 21 + w. The advantages of this 
choice will become clear in what follows. In the new coordinates (J, t) lemma [l] rewrites: 

Lemma |l|' For a fixed Lq ^ and given an initial condition (Jo,£o) with Jo > 0, its iterated by the 
billiard map, (Ji,ti), is given by the formulae: 

Ji = Jo-(w + w 1 ) + 2 R }~f° (6) 

ii — to 

ALl + {J -w f R\-I% 

{tl to) arJ Jo w ° + 17^77' (7) 

Since R(t) is at least C 2 , there exist £o and £i such that 

R\ - «l = - 1„) + «,)' - u ^(t„ - tl f, (8) 

For Lq 7^ and by inequality (||) we have 



?7? 2 

|Ji-Jo|<W— (9) 



where w' = max{\w'(t)\,t € [0,T]}. 



Remark: Suppose that the breathing circle moves slowly, i.e., R 2 (t) = 1 + eg(t), for a small parameter 
e and g a C k , T-periodic function. Then w'(t) = eg" it) — 0(e). Let J n be the value of J at the n th 
impact. For n = O(-) 

or? 

\J n - M < nw'-— = 0(1) 
which shows that J is an adiabatic invariant (see |l| for the definition of adiabatic invariant). 
Let B x denotes the cylinder (A, +oo) xS ! = {(J, t),J>\,te S 1 }. 

Lemma 3 Suppose that R(t) is a strictly positive T-periodic C k function, k > 2. Let 

M : B -2 — > M (B ,-=2 ) C B^., 

w "'T^oT ; \ [Zol ^ > 

be defined by formulae ^) and (Qj. Then M is a C k ~ 2 diffeomorphism, preserving the measure 
J 

AH 2 



dfl = -drzdJdt 



Proof: For a fixed Lq ^ 0, we have automatically from (^) that 

27? 

J > w'(— — r + 1) =>• Ji > w' > 
\ L o\ 

The proof follows immediately from this fact and lemma H. 
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Lemma 4 Given any curve T in the cylinder B_ —2 homotopic to S 1 x {0} then M(T) (IT ^ 0. In 
other words, M has the property of intersection. 



Proof: M preserves a measure, absolutely continuos with respect to the Lebesgue measure. This implies 
the property of intersection . ■ 



4 Invariant Curves 

The purpose of this section is to investigate the behaviour of the breathing circle billiard map in the 
neighbourhood of 00. We shall show that a map can be defined which exhibits some properties of twist 
maps. More specifically, using a result proved by R. Douady [Q we will prove the existence of invariant 
curves for values of J sufficiently large. 

For a fixed L$ ^ 0, and since 

2if 

J > w'{—— + 1) => Ji > v? > 
\ L o\ 

the change of variables (z,t) = (j,t) is well defined on B_ r 2 -2 and is C°° . 



Let < <5 < — 1^ q| . Then the map 



w'(2R 2 + \L a \)' 



N : (0,5) xS 1 M+xS 1 
(z ,t ) h-» M(l/z ,t Q ) 



defined by 



1 R 2 — R 2 

z i = z ° 1 , u with 9(ti,t ) = -(w + 101) + 2-i (10) 

l+g(ti,t )zo h~t 

[AL 2 z 2 + (1 - w a z ) 2 ] (h - t Q ) + m 2 z (w z a - 1) = 4R 2 z 2 R }~^ (11) 

t\ — Iq 

is a C fc_2 diffeomorphism, with the property of intersection. 

Lemma 5 N can be extended to a C k ~ 2 diffeomorphism N on a neighbourhood of {0} x S 1 . N has the 
property of intersection and a development of the form: 

zi = z + 0(z 2 ) 

h= t + AR 2 zo + 0(z 2 ) 

Proof: First of all, let us show that t\ = ti(to, zo) can be extended for zq < 0. R(t) being C k , k > 2 then 
f{h — to) — R\ - R 2 - w {ti - t ) is a C fe_1 function with /(0) = /'(0) = 0. By Hadamard's Lemma 
f(t\ — to) = [ti — to)h(t\ — to), where ft, is a C k ~ 2 function. Using ( |TT| ) and discarding the trivial solution 
t\ = to, we have: 

[4L 2 z 2 + (1 - w z ) 2 ] (h - t ) + 4Rlz (w zo - 1) = 4R 2 z^ [w + h(h - to)] . 
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Let 

F(t , h,z ) = [ALlzl + (1 - w z ) 2 ] (ti - t ) + mlz (w z a - 1) - ARlzl [w + h(h - t )\ . 

F is a C k ~ 2 function on S 1 x S 1 x R such that F(t , t , 0) = and §^(t , t , 0) = 1. If fc > 2 and for each 
fixed to G S 1 , by the Implicit Function Theorem we have t\ = t\(zo,to) on a neighbourhood of (0,to)- 
Moreover, as S 1 is compact and ti(0, to) = £o for every to, we can find a neighbourhood of {0} x S 1 where 
t\ = ti(zo,to) is a C k ~ 2 function. 

This implies that g(z ,t Q ) = -(w(t ) + w(ti(z ,t ))) + has a C k ~ 2 extension on a 

neighbourhood of {0} x S 1 . Clearly it is bounded and, perhaps on a smaller neighbourhood, we can take 
\g(ti,to)zo\ < 1. So, Zi(zoi^o) has a C fe_2 extension on a neighbourhood of {0} x S 1 . 

Analogously, A^ 1 can be extended on a neighbourhood of {0} x S 1 . N is then extended to a C k ~ 2 
diffeomorphism TV on a neighbourhood of {0} x S 1 . 

Since the formulae remain unchanged, the property of intersection is preserved for curves on zq > or 
z < 0. But N(0, t ) = (0, to) for any t e S 1 . Then, if a curve T cuts the circle {0} x S 1 , N(T) n T ^ 0, 
and TV has the property of intersection. 

Finally, a simple calculation leads to the development of N. m 

Proposition 6 If k > 7, we have a family of invariant curves for N that protect {0} x 5 1 , in the 
following sense: Every orbit of N with no intersection with {0} x S , stay at a bounded distance, above 
and bellow, from {0} x S 1 . Those bounds are arbitrarily small strictly positive constants. 

Proof: Lemma H give the conditions to apply a corollary of Herman's Theoreme des Combes Translatees, 
proved by R. Douady ( ||, page III - 8) and the result follows. ■ 

5 Boundedness of the velocity 

Coming back to the coordinates (J,t) and the diffeomorphism M, let us denote by (J n ,t n ) the points 
of the orbit of (Jo, to) under M. An admissible initial condition for M is the set of (Jo, to) such that 
J„ > 0,Vn (which is the condition of existence of the next impact given in lemma [l]'). 

Theorem 1 Given the billiard on the breathing circle x 2 + y 2 = R(t) 2 , with R(t) a strictly positive T- 
periodic C k function, k > 7, then, for any admissible initial condition, a particle will move with bounded 
velocity. 

Proof: For Lo ^ 0, since N has a family of invariant curves that protect {0} x S 1 and N = N for 
zo > 0, we have a family of invariant curves of M that protect oo, i.e., there exist a family of functions 
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r„ : S 1 — > ( w '(jj^ + approaching infinity with norm C k 2 , whose graphs are invariant under 

M. 

Given {Jo, to) such that Ti(t ) < J < T m (to) for some / and m, we will have that < mi teS iTi(t) < 
J n < snp tGS iT m (t), for every n. Therefore (Jo,io) is an admissible initial condition and J n will remain 
bounded. 

On the other hand, if (Jo, to) is bellow all the invariant curves, then J„ will be bounded above although 
it may be equal to zero. Anyway, for every admissible initial condition, J„ will remain bounded. 

Since J n = 2(I n + ^f-) = 2(—R n v n ■ fj n + it is obvious that J n is bounded if and only if v n — \v n \ is 
bounded. 

If Lq = 0, we have the motion along the diameter of the breathing circle, that corresponds to the 1- 
dimensional Fermi accelerator, which has bounded velocities for any admissible initial condition (see, for 
instance, Q]). ■ 

Remark: This result is also true if we consider the billiard in any higher dimension breathing sphere 
5", instead of the circle S , since the conservation of the angular momentum implies that the movement 
will occur on a plane and so, the effective boundary will be the breathing circle. 

In terms of the dynamics of the breathing circle billiard, we can translate Theorem 1 in the following 
way: For each fixed Lq, there exists a bound K, depending on R(t) and decreasing with Lq, such that for 
J > K, there are infinitely many spanning curves that may not exist for J < K. This fact is illustrated 
on the pictures bellow where we plot the phase-space for R(t) = 1 + jjp cos 27rf and different choices of 
Lq. (Because of the special form of R(t), J is, in fact, an eternal adiabatic invariant and, for big J's, all 
curves are invariant.) 




Figure 1: Phase-space of the breathing circle for t 6 [0, 1), J € [0, 8], and Lq — 0.25, 0.98 and 1.3. 



Clearly our approach does not apply to the low energy regime (i.e., for small Lo' s and/or small J's). 
For the 1-dimensional case, a simplified model was introduced by Lichtemberg and Liebermann Q. This 
approximation and some comparisons with the standard map seems to give some valuable insights about 
this energy regime. In the simplified model one allows the boundary to interact with the particle through 
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momentum exchange (via a given w(t) periodic in time), but assuming that the boundary does not change 
in time and so R(t) = Rq- In this case, the associated billiard map is given by 

J% = Jq- (w + Wl) 



4 -Rn 

h = t + 2 . ^(Jo - wo). 

4Lg + (J - w y 

This simplified model allows unphysical negative values of J, which are forbidden in the complete model. 

On the next set of pictures we show the phase-space of the simplified model with the same conditions 
and w(t) we took in figure Except for a small neighbourhood of J = 0, the dynamical behaviour seems 
astonishingly similar. 




35«S*!».«. » ,-, ZZ.TZ Z ... i»J(E5»55 



Figure 2: Phase-space of the simplified breathing circle for t G [0,1), J G [0,8], and Lq = 0.25,0.98 and 
1.3. 

Finally we point out that the existence of a conserved quantity (the angular momentum) is the clue of 
our proof, since it reduces the dimension to two. Unfortunately, it is easy to prove that a 2-dimensional 
moving billiard, with convex frozen boundary and with only normal deformation, preserves the angular 
momentum if and only if it is the breathing circle. So, the extension of this result to other time dependent 
billiards will depend on the existence of some constant of motion or to the extension of the techniques to 
higher dimensions. 
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